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I I review some recent work on the problem of multiparticle production in a 0^- 

QhI theory with an 0{N) symmetry. Threshold amplitudes with fixed number of pro- 

Q^' duced particles are exactly calculated at large- to all loops and vanish on mass 

^ ■ shell for 2^n when n>2 due to a dynamical symmetry. I consider an extension to 

the cases when the 0{N) symmetry is softly broken by masses including sponta- 
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neous breaking of a remaining reflection symmetry. The exact solutions are obtained 
by the Gelfand-Dikii technique of finding the diagonal resolvent of the Schrodinger 
5^ I operator which emerges due to factorization at large N. I report also some new 

results on the diagonal resolvent for a general Poschl-Teller potential, which could 
be useful for calculations of multiparticle amplitudes in the standard model. 
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1 Introduction and resume 



1.1 n\ 

The history of the subject began in 1990 when Ringwald [jl| and Espinosa utihzed 
n\, which is due to n identical scalar particles in a final state, to get a large amplitude 
of B+L violating processes in the standard model due to instantons in the dilute gas 
approximation. Soon after Cornwall and Goldberg applied a similar idea to a 
problem of perturbative multiparticle production in the 0^ theory with a small coupling 
constant A. Some other early references include also ||, |^. 

In the 0^ theory one considers the process when a virtual particle with the 4-momen- 
tum {nm, 0) produces n on-mass-shell particles of mass m at rest with the 4-momentum 
(m, 0) each. Such a process may appear in a high energy collision of quark and antiquark 
when n Higgs particles are produced. Since the amplitude of the process in proportional 
to n\: 

a{n) oc n\ , (1.1) 
a naive estimate of the cross-section in the tree approximation is 

a,,-^nH~^|a(n)|2~n!A"-S (1.2) 

which may become large for n > 1/A. 

A special comment is needed about "triviality" of the 0^ theory. If the 0^ theory were 
a fundamental one, the renormalized value of A would vanish. A modern look at this 
problem is that it is no longer fundamental at the distances 

r~m-^-^^'/=^^ (1.3) 

which are very small if (renormalized) A is small and play a role of a cutoff. Therefore, 
the estimate ( |1.2| ) makes sense for small enough A when n <C exp (Stt^/SA) can still be 
large. 



1.2 The tree level 

The explicit results for the tree amplitudes at threshold in the 0^ theory 

/ . \(n-l)/2 

a(n) = -{n' -l)n\m'\—-\ , (1.4) 

where c = 8 |Q for an unbroken refiection symmetry (j) —>■ —0 and c = 2 for a broken 
one, demonstrate the expected factorial growth. These results are obtained solving exact 
tree-level recurrence relations between the amplitudes a{n) with different n, which is an 
alternative to semiclassical evaluations at large n. A bridge between the two approaches 
relies on a functional technique |]lOl which turns out to be very useful to go beyond the 
tree level. 
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1.3 Unitarity violation 



While the threshold amplitudes have a vanishing phase volume, they are a lower bound for 
tree amplitudes with nonvanishing spatial momenta of produced particles This bound 
looks like the r.h.s. of Eq. (|1.4| ) with m substituted by the variable uj = E — {n — l)m 
where E is the total energy of produced particles in the rest frame {uj = mii all particles 
are produced at rest) and c = 9orc = 8/3 [jlT], |l^ for the unbroken or broken symmetry, 
respectively. 

This behavior of the tree amplitudes leads to a violation of unitarity |P, 11, 12, T3| . 



Moreover, if one reconstructs the perturbative expansion by the unitarity relation, it 
means under some natural assumptions that the perturbation theory breaks down starting 
from some finite order. Therefore, loops should be taken into account. 

1.4 The nullification 



The one-loop calculation of the l^n amplitude by Voloshin [|14| discovered a very in- 
teresting property of the tree 2— threshold amplitude which vanishes for n>4 when 
incoming particles are on mass shell. This nullification occurs for dynamical reasons be- 
cause of the cancellation between diagrams. It has been extended [15-25] to more general 
models and is crucial for calculations of the amplitude 1— i>n at the one-loop level. A 
dynamical symmetry which may be responsible for the nullification has been discussed in 
Ref. EH]. An interesting question is which properties of the tree and one- loop amplitudes 



could survive in the full theory which includes all loop diagrams. 
1.5 0(oo) symmetry 

The simplest exactly solvable to all loops model is the 0(A^)-symmetric 0^ theory in the 
large- limit when only bubble diagrams contribute if there is no multiparticle production. 
The amplitudes of multiparticle production at threshold are calculated for this model in 
the large- A^ limit at fixed number of produced particles n in Ref. . The results, which 
are reviewed in Sect. 2 below, are quite similar to the tree amplitudes while the effect of 
loops is the renormalization of the coupling constant and mass. The form of the solution 
is due to the fact that the exact amplitude of the process 2^n vanishes on mass shell 
for n>2 at large A^ when averaged over the 0{N) indices of incoming particles. An 
extension to the cases when the 0{N) symmetry is softly broken by masses including 



spontaneous breaking of a remaining reflection symmetry is considered in Sect. 3. 
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1.6 The rescattering problem 



Since the large- iV multiparticle amplitudes are similar to the tree-level ones, they are not 
very instructive from the viewpoint of unitarity restoration by loops (while agree with 
unitarity by themselves). More essential for this are loops associated with rescatterings of 
the produced particles which are controlled at large n by the parameter [|l^, 0. As 
A ~ they give corrections to the large- amplitudes which are controlled by /N. 

The rescattering contribution is expected to be calculable by semiclassical meth- 
ods [18, 29|. Nice explicit results showing how unitarity is restored by this mechanism are 
obtained fS^] in 2 + 1 dimensions by an alternative method of summing infrared logarithms 
with the aid of renormalization group. The conjecture about exponentiation of the rt^X 
(or /N) corrections, which is confirmed by various considerations [31], can help to solve 
the problem. 



1.7 Yet other app Heat ions? 

The techniques developed in studying multiparticle production might be useful for other 
problems (some of them are mentioned in Sect. 4). An idea is that for the given kinematics, 
which reduces the system at the tree and one-loop levels or at large A^ to an integrable 
lower-dimensional one, calculations could be simpler than, say, for all momenta to be of 
the same order. 

Equations of this type appear in many problems with an emission of energy. In 
particular, the large-A^ Schwinger-Dyson equations considered in Subsect. 2.2 coincide 
with those describing a non-equilibrium time evolution of quantum systems.^] The 



difference is that the latter problem is in real rather than in imaginary time and, what is 
more essential, the boundary conditions are quite different. 

2 Unbroken O(cxd) symmetry [27] 
2.1 The definitions 

The 0(A^)-symmetric 0^ theory is defined in the 4-dimensional Minkowskian space by the 
Lagrangian 

c = 1(9,0") (9,0^) - |m2(0y) - \H<pyr , (2.1) 

where b = 1, . . . ,N and the summation over repeated indices is implied. 

The amplitudes a|j^...j,^(n) of production of n on-mass-shell particles with the O(A^) 
indices 6i . . . 6„ at rest by a (virtual) particle with the 0{N) index b and the energy nm 

I thank E. Mottola for this comment. 
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are given by the generating function P, ITO 



n>l 



ra! n"' — 1 



n>3 

where r = ixo is the imaginary-time variable and 

a\n) = al..,^in)e'...e-. (2.3) 

The multiphcation by the 0{N) vectors C,^' results in a symmetrization over the 0{N) 
indices of the produced particles. The tree-level estimate gives 

a\n)^{xe) ' e. (2.4) 

Since A ~ in the large- limit, we choose ~ for all the amplitudes ( p.4|) to be 
of the same order in 

It is convenient to introduce the coherent state 



E^e'^-'-e'" (&l---&n| . (2.5) 



Then Eq. (|2.2| ) can be written as 

^\T) = (E\<l)\6,-tT)\0) . (2.6) 



The second quantity which appears in the set of the Schwinger-Dyson equations of 
the next subsection is the amplitude D^^,,,f^^{n;p) of the process when two particles a and 
b with the 4-momenta p + nq and —p, respectively, produce n on-mass-shell particles with 
the 0{N) indices 6i . . . 6„ and the equal 4-momenta q = (m, 0) in the rest frame. We 
define again 

D^\n;p) = DtU^in;p)e'...e- (2.7) 
and introduce the generating function 

GZ\r,r') = /^e^^^-^(S|0'^(f,-^r)0^(O,-zr')|o)-<5(3)(p-)$"(r)$^(r') 



i2n) 

) 

2u 



27r „ n 



n=0 



Here 



CO = \Jp + w? , (2.9) 

p = {e,p) and the subtraction cancels disconnected parts which are not included in the 
definition of D°'^{n]p). The meaning of C^ij^r') is that of the Green function in the 
imaginary-time representation while D"'^{n]p) is the one in the energy representation. 
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2.2 Scwinger— Dyson equations 

The Schwinger-Dyson equations can be extracted from the following identity 



6S[4>] 



50" (a;) 



m 







6F[<P] 







(2.10) 



where 5* is the action of the model and F[(j)] is an arbitrary functional of (p. Eq. (|2.10|) 
results from the invariance of the measure in the path integral under an arbitrary variation 
of Substituting F = 1, one gets from Eq. ( p.lOj ) 



+ 



- I -ir)|0) - A -zr)0'^(f, -ir)|o) = . (2.11) 



In the large- limit one can split the matrix elements of the operator 



as 



bi ...bn\ 



/,2 ial 



o) = E E 



,p' ni+n2=n ^l-^2- 



p(6i) . . .p(6„J|0^|O (p(6„,+i) • • •p(fen)iriO) 



(2.12) 

where p and p' stand for permutations. This formula holds for n finite as — * cxd and 
extends the standard factorization of 0(A^)-singlet operators at large to the case of 
multiparticle production. According to the definition ( p.5| ), we rewrite Eq. ( p.l2| ) in the 
form 

'E\<f)\x, -zT)r{x, -tT)\0) = (S|02(f, -,r)|0) {E\r{x, -^r)|0) . (2.13) 



Using Eqs. (|2.6|) , ( 2.13 ) and translational invariance, Eq. ( P . 1 1 ) can be written as 
where we have denoted 



^;(r) = A S|0^(O,-zr)|O) = A$^(r) + 



A 
2^ 



duo uVlo^ - m2 CJiT, t) (2.15) 



and the last equality holds due to the definition (|2.8| ). 

The meaning of the amplitude G"" which is summed over the 0{N) indices of incoming 
particles becomes clear after the obvious decomposition 



acb 



G'^(r,r') + ^Gf(r,r'). 



:2.16) 



At large A^ one gets 



(r, r') = (l - ^) Glir, r') + ^Gf (r, r') 



(2.17) 



Only enters the Schwinger-Dyson equations at large A^. 
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In order to close the set of equations, one derives an equation for the propagator 
G'^{t,t') which corresponds to F = 0°(O, —ir') in Eq. ( p.lO| ). Using the large-A^ factor- 
ization and integrating over d^x, one gets 



) = -Sir - r') 



(2.18) 



The set of Eqs. ( |2.14|) , (|2.15|) and ( |2.18|) is closed and unambiguously determined the 
solution presented in the next subsection. 



2.3 The exact solution 



Eqs. (|2.14| ), (|2.15|) and (|2.18|) look like a quantum mechanical problem which can be 
solved applying the Gelfand-Dikii technique which is described in Subsect. ^A. 
The solution, which satisfies proper boundary conditions as r —>■ —(yo, reads 



vr{t 



2u Auioj'^ — m^) ' 



where 



and 



A 



R 



XnN ' 



<l>"(r) 



I _ MIq 2miiT 



(2.19) 

(2.20) 
(2.21) 

(2.22) 



The solution is expressed via the renormalized mass and coupling which are 
related to the bare ones m and A by 



m 



m 



R 



XN 

47r2 Jm\ 



duj \Juj'^ — m 



and 



1 _ 1 
A ~ A^ 



uo^ — ra 



R 



^2.23) 



(2.24) 



which exactly coincide with the standard renormalizations of mass and the coupling con- 
stant at large A^. The renormalization of v{t) is given by 



'Vr{'t) = v{t) + — m 



R ■ 



(2.25) 



Eqs. ( ^.19| ), ( |2.22| ) coincide with the large- limit of the tree- level results for G^{t, t) 



|[T7| and $"(r) ||TOl in 0(A^)-symmetric 0^ theory, while the only effect of loops, besides 
the renormalization, is the change of Xr given by Eq. ( p. 211) . 
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a) b) c) 



Figure 1: The tree diagrams for D-^(4;p) at large N. The continuous hnes are associated 
with propagation of 0{N) indices. The contributions of the diagrams a), b) and 
c) equal, respectively, —1/4, 1/8 and 1/8 for the given on-mass-shell kinematics 
when the momenta of the incoming particles (in the units of mass) are (2, \/3, 0, 0) 
and (2, — \/3, 0, 0), while that of each produced particle is (1, 0, 0, 0). The sum of 
three diagrams vanishes which illustrates the nullification of 2^4 on mass shell 
at the tree level. 

The fact that (|2.19| ) has a pole only at uj'^ = m\ means the nullification of the on- 
mass-shell amplitudes D^{n\p) for n>2. An analogous property holds |^ at = 1 



for the tree level amplitudes in the case of the sine-Gordon potential and for the (j)^- 
theory with spontaneously broken symmetry It can be explicitly demonstrated by 
the cancellation of the tree diagrams for the process 2— which are depicted in Fig. |^. 

The large- amplitude ( |2.22| ) is real which is due to the nullification. This is different 
from the A^ = 1 one-loop result |]l4l while is similar to the case of spontaneously broken 
reflection symmetry where the 1— *-n amplitude is real to all orders [Q. The fact that the 
amplitude 1— is real at the one-loop level is illustrated by Fig. |^. While the sum of the 
imaginary parts vanishes, the real part is described by Eqs. ( |2.21|) , ( |2.22| ). 



Analogously, the imaginary parts of the one-loop diagrams for the 2— i>4 on-mass-shell 
amplitude -D^(4; p) are mutually cancelled for the same reason while the real part is given 
by Eqs. ( p.l9|) , ( p.21| ) and ( |2.22|) . This amplitude is also real to all loops. Therefore, both 
amplitudes are real due to the nullification which, in turn, is due to a hidden dynamical 
symmetry. 



3 Softly broken 0(oo) symmetry [28] 
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a) b) c) 

Figure 2: The one-loop diagrams for a^{5;p) at large A^. The dashed lines represent unitary 
cuts over two-particle intermediate states. The imaginary part vanishes after 
summing over all three diagrams since the amplitude D^(4;p) vanish on-mass- 
shell at the tree level as is illustrated by Fig. |l|. 

3.1 The 0{Ni)xO{N2) model 

The model consists of an 0{Ni) multiplet x''; b = l,...,Ni, coupled to a 0(A^2) scalar 
field 0^, P = 1, N2. The Lagrangian reads 

£ = l(d,x')id,x') + lid,<P^){d,<P^)-krr^l{x\'^ (3-1) 

where interaction term is 0{Ni+N2) invariant while the mass term possesses the 0{Ni) x 
0{N2) symmetry for nonequal masses mi and 1712- All the formulas of this section are 
proper extensions of those of Sect. 2 to mi 7^ 1712. 

We denote by ^^^^^^ ^^(n, k) the amplitude of production of n on-mass-shell parti- 
cles X with the 0(A^i)-indices 61 . . . 6„ and of k on-mass-shell particles (f) with the 0{N2)- 
indices Pi, . . . , Pk at rest by a (virtual) particle x with the 0(A''i)-index b and the energy 
nmi + km,2. We also denote by B^_^ ^^{n, k) the amplitude of production of n on- 

mass-shell particles x with the 0(iVi)-indices 61 . . . 6„ and of k on-mass-shell particles 
with 0(iV2)-indices /^i, . . . , at rest by a (virtual) particle with the 0(iV2)-index P and 
the energy nnii + km2. Analogues of Eq. (|2.3| ) are 

A\n, k) = (n, k) ei' . . . ■ ■ ■ , (3.2) 

B^{n, k) = K....„,ft...A {n, k) ei' ■ ■ ■ ■ ■ ■ il' . (3.3) 



An analog of the coherent state ( p.5|) is 

00 1 

(S| = H {hi... hn,Pi...Pu\ , (3.4) 

n,k=0 
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while the generating functions for the amphtudes (p.2|) , (|373|) are 

vl/"(r) = (s|0"(O,-zr)|o)= 7712^2 e'"-+ J] ' ^, , 2Y^>2- (3-6) 



n,k>3 

The definition of the generating functions C^i^jj-.T^r') and C^^ {ljJ]t,t') which extend 



n\k\{{nmi + km2Y ~ ^2 

and (|2.8| ) is obvious. 

The Schwinger-Dyson equations which extend Eqs. (|2.14|) , (|2.15| ) and ( p.l8|) read 

|^-m?-7;(r)|$«(r) = 0, (3.7) 

|^-m^-t;(r)|vl/"(r) = 0, (3.8) 

- a;^ - 7;(r)| Gj(c.; r, r') = -5(r - r') , (3.9) 

1^ - a;^ - 7;(r)| GJ(c.; r, r') = -5(r - r') , (3.10) 



where 



v{t) = X (h|(x'(0, -zt) + 02(0, -zr))|0) = A($2(r) + vl/2(r)) 



+ ^^y duJUJ^/uJ^^^^Gl{uj;T,T) + ^-^ J duu^u^ - Gj(^;r,r). (3.11) 
The derivation uses factorization at large A^. 



3.2 Unbroken reflection symmetry 

The solution to Eqs. (|3.7|) - (|3.11|) is very similar to that of Subsect. 2.3, while there are 
important differences. It reads 



^T. ^_ 1 Ai$2 X2<f' 



where 

= Ai$2 + A2^^ (3.13) 
Ai = , A2 = , (3.14) 

and 



2 

m2H 



(3.15) 
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Figure 3: The one- loop diagram for the process x + X 
part for m2R > miR given by Eq. ( ^.21 ). 



which possesses an imaginary 



2 

IR 



m 



2R 



(3.16) 



while and are similar to those of Ref. \24 



K 



m 



2R 



/3 _ 



l + 2Ar 



K 



m 



2 ^1 

IR / 



2Ai ^2 


2X2 


2 ^1 




T^lR 


^Ir 


2Ai 2 


2X2 


2 ^1 




^IR 


^Ir 



Zo + A1A2 



^ 2 2 

2 2 ^1 ^2 
^lil^2_R / 



+ A1A2- 



2 2 
"^1^2 



1R"''2R. 



with 



and 



K 



mm - m2R 
niiR + m2R 



(3.17) 
(3.18) 

(3.19) 
(3.20) 



The constants Ai and A2 are complex in general. Assuming for definiteness that 
> one can easily see that the constant Ai is real while the imaginary part of A2 
is given by 

I,„i = -^^5H^. (3.21) 



A2 Stt 2m2i? 

The appearance of the imaginary part is because the masses of the 0{Ni) and 0(A^2) 
particles are not equal. For mlj^ > mm the imaginary part is associated with an inelastic 
process X + X + which is given to lowest order by the diagram depicted Fig. p. 



3.3 Broken reflection symmetry 

For N2 = I when is a singlet field, the symmetry of the Lagrangian ( p.l|) under the 
reflection (p — > —(f) is spontaneously broken if < 0. The remaining 0{Ni) = 0{N~1) 
symmetry is unbroken for ml > — I^TiH (for mf = m\ < the whole 0{N) symmetry were 
spontaneously broken down to 0{N—1) and the associated massless Goldstone bosons 
would appear). 
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The vacuum expectation value of the field and the physical mass of the field x ^'^^ 

\m2R\ 



^2 5 



while the physical mass of the shifted field 0' 



(3.22) 

?7r is determined self-consistently by 

(3.23) 



with Ai, A2 given by Eq. ( p.l4| ) and 
h = N 
I2 = N I dujJuj 

.1 m • 



(cij^ — 771^/4) (cu^ — m^) uo'^ 

1 1 
oj"^ — m?/4 oj"^ 



(3.24) 
(3.25) 



The parameters Ai and A2 are real in contrast to the case without the spontaneous breaking 
of the refiection symmetry. 



The solution reads 



G^(^;r,r) 



00 



Aicj$^ 



2(cu2 - my A) 4(cu2 - m2/4)(cu2 - m2) Auj{uj'^ - ml/A) 



t;^ = Ai$^ + A2^^ - -r^ 



2 



where 



^ _ (2mx-m.0) ^2 
(2mx+m,0) Irin 



Y £2 8A1 ^2 _|_ 4Ai (2m.x-m^) ^ ^2 

2:2 



I6A1 



= Vr + Vr 



1 ~^ 



1 - 7^ 



8A1 ^2 1 4Ai (2m^-m^) ^ 2 



and 



4 



2:2 = rn^^2e' 



(3.26) 
(3.27) 

(3.28) 
(3.29) 
(3.30) 



This solution is similar to the tree level one [241 and is real. 



The poles of (|3.26|) at u = m^ and u = m(^/2 are associated with production of two 
X-particles or a single particle, respectively. There is no pole at = which signals 
that the amplitude for the production of two particles vanishes in contrast to the cases 
without the spontaneous breaking and A'^ = 1 with broken refiection symmetry [|1^]. It is 
easy to see shifting \1/ = 77/? + that ( p.26| ) and ( p. 271 ) recover the free values l/2uj and 
0, respectively, to zeroth order of perturbation theory in A^j around the shifted vacuum. 
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4 Generalizations 



4.1 Gelfand— Dikii technique for Poschl— Teller potential 

The results presented in Sects. 2, 3 are obtained by the Gelfand-Dikii technique which 



represents the diagonal resolvent of the Schrodinger operator in Eq. (|2.18 

Ri[v] 



as 



1=0 



21+1 



where the differential polynomials Ri[v\ are determined recurrently by 



Ri+i[v] -- 

and the inverse operator is 



1 ( i n2 

2 12 



- V - D-^vD) Ri[v 



Rn 



(4.1) 

(4.2) 
(4.3) 



D-'f{T) = / dx fix) . 

J — oo 

Eq. ( [4.2|) stems from the fact that R^^ [v] obeys the third order linear differential equation 

^ ^]D'^-Dv- vD) R^ [v] = uj'^DR^ [v] . (4.4) 



2 V2 



The potential ( p.20|) , (|2.22| ) is a particular case of a more general Poschl-Teller poten- 



tial 



v{t) = 



sis + l] 



k{k + l) 



sinh^ (m(r — To)) cosh^ (m(r — Tq)) 
with s = 1 and k = 0. The simple answer ( ^^.191 ) for the diagonal resolvent is because 

1 



(4.5) 



Ri[v] = m^'-^Ri[v] = --m^'-^v for s = l,k = 0, 



(4.6) 



i.e. Ri[v] with / > 1 reduces to -Ri[f]. 

It is easy to extend Eqs. ( p.l9| ), ( p.20| ) to case s>l,/i; = 0(or/i;>l,s = 0) when 



RM = E 



Luiu"^ — rn?) ■ ■ ■ (u;^ — Prn?) 
and the coefficients Ci are determined recurrently by 

/(/ + 1) -s(s + l) 



i+i 



l + \)Ci. 



(/ + l)s(s + l) 

These formulas can easily be proven substituting into Eq. ( [4.4| ) and noticing that 

/(/ + 1)-S(S + 1) ^, , 



1 n n2 



2 V2 



-V- D-^vD ) v' = t^m'v' + 



.1 ?2_2„; 



(4.7) 



(4.^ 



(4.9) 



It is evident from Eq. ( [4.9|) that f*, which is the highest power of v in (|4.7| ), is an 
eigenvector of the operator entering Eq. ( [4. 21) and, therefore, -Rz[f] for / > s reduces to 
RAv\....Rs\v\. 
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The simple form ( [4.7|) of the diagonal resolvent for the potential ( |4.5| ) with = is 



due to the fact that the spectrum is very simple: there are s poles at = 1, . . . , s (in the 
units of m). For nonvanishing both s and the spectrum is slightly more complicated: 
there are M = max{s, k} poles which go from 1 up to \s — k\ with the step 1 and then 
from |s — A;| + 2uptos + /c with the step 2. Say, the spectrum for s = 2 and k = 5 
is 1,2,3,5,7. This spectrum can be obtained from the one for = (or s = 0) by a 
supersymmetry transformation^] which relates {s,k) with [s + l,k + l). It was explicitly 



done for k = 1 and any s by Smith |25|. 



For the potential ( [4 .51) Ri[v] with I > M = max{s, /c} reduces to Ri[v], . . . , Rm[v\. 
This is based on the fact that 

\ (^D^ -V- D-^vD) sinh-2' r cosh'^'^ r = (s + kf smhr"^' t cosh"^^ r (4.10) 

for the potential ( [4.5| ) (I put m = 1 and tq = for brevity), i.e. there exists an eigenvector 
of the operator entering the recurrence relation ( |4.2| ). As is shown in Ref. [^], if $(r) is 



an eigenvector of the Schrodinger operator on the l.h.s. of Eq. ( |2.14| ) with an eigenvalue 



E and $(— oo) = 0, then $^(t) will be an eigenvector of the operator on the l.h.s. of 
Eq. ( ^lOD with the same eigenvalue. Therefore, spectra of the two operators coincide. 



Since the spectrum is known, a modification of ( [4.7|) for k ^ is obvious while 
awkward. An analog of Eq. ( [4.6|) can be obtained as follows. For the given spectrum 
mi, . . . , mjvf let us calculate 



M AI 



V^'^^\u) = n(^' - ^- ) = - E -'^cu^^*'-*) (4.11) 

i=l i=l 

to determine the coefficients d\^'^\ Then, the reduction of -Rj\/+i[i;] reads 

M 

Rm+1 = Dj' Rhi-i+i- (4-12) 



The reduction of Ri[v\ for / > M+1 can be obtained from Eq. ( |4.12|) recurrently using 
Eq. (13). 

Let us illustrate Eq. ( |4.12| ) which relies on Eq. ( [4.10| ) by an example of (s=2, k=b) 
with the spectrum mentioned above. One gets 

r'(2.5)(^) = (^2 _ ^^(^2 _ 4) (^2 _ g)^^2 _ (^2 _ ^ 

u^^ - muj^ + 2310cj6 - 20812u;^ + 62689u;^ - 44100 (4.13) 



and 



i?6 = 88i?5 - 2310/?4 + 20812i?3 - 62689i?2 + 44100/?i , (4.14) 



§ The potential (15) with k — \ and arbitrary s (which is related to the ratio of fermion to vector-boson 
masses) appears in the standard model for longitudinally polarized vector bosons ||2^, . 
^ I am grateful to E. Argyres and A. Johansen for discussions of this point. 
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a) b) c) 



Figure 4: Some large- diagrams for the symmetrized amplitude G^{4;p). 



while the diagonal resolvent is 

R U,! — J_ _|_ fii _|_ R2-R1 _|_ fi3-5R2+4fli 



I -R'.4~14fl3+49-R2-36Ri i /;s-39/;4+399-R3-1261R2+900-Ri (A^'^\ 

"^c<;(l^2_i)(^2_4)(^2_9)(^^2_25) "T <^(a;2_i)(a;2_4)(aj2_9)((^2_25)(^2_49) • l^i--^^; 

The coefficients of the numerators in the last term on the r.h.s. are just -Dp'^"* associ- 
ated with the spectrum 1,2,3,5 for (s=l, fc=4) as is prescribed by the supersymmetry. 
Analogously, the numerators of the 5th, 4th and 3rd terms on the r.h.s. correspond to 
the (s=0, fc=3), (s=0, A;=2) and {s=0,k=l) spectra, respectively. Substituting (^4.5|) and 
using Mathematica, ( [4.15|) can be written explicitly as 



p r 1 _ J_ I 3sinh~2T(4+cosh~2T) , 32 sinh~"^ t(21-55 cosh~2 r+35 cosh~^ r) 
— 2cj + 2uj{ui^-l) + 2aj{w2_i)(^2_4) 

. 32.52 sinh''^ T cosh'2 r(10-35 cosh~2 ,-+28 cosh-"^ r) 
+ 2a.)(a;2-l)(w2_4)(^2_9) 

I 33-52-72 sinh~^rcosh~^T , 3" -52 -72 sinh~"^ r cosh~^° t (AAR\ 

"^2a;(Lj2_i)(^2_4)(^2_9)(^2_25) "T 2c<;(w2_i)(^2_4)(^2_9)(j^2_25)(a;2_49) • l^^.iU; 

The highest vector is in agreement with Eq. ( [4.10| ). 



4.2 Symmetrized amplitude and renormalons 

The symmetrized amplitude 

D'in;p)=D'^\n;p)^ (4.17) 

reveals, on the contrary, a surprisingly nontrivial behavior even at large N. Some of the 
diagrams which contribute to D'^{4;p) are depicted in Fig. ^ The combinatorics is now 
different from that of the diagrams of Fig. |l] so there is no cancellation of tree diagrams 
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a) b) c 



Figure 5: The diagrammatic representation of some diagrams for D^{A;p) at large N inte- 
grated over d^p. The diagrams b) and c) are of the renormalon type. 



for D'^{A;p) on mass shell which happens for n>6 according to the explicit result [17| for 
D^{4:;p) at the tree level. Each of the vertices can be "dressed" with bubble chains and 
each of the lines of outgoing particles can be substituted by the exact amplitude a''(ra) to 
obtain a diagram with more produced particles. 

One can integrate D^{4:;p) over the 4-momentum p to obtain diagrams of the vertex 
type. Some of them are depicted in Fig. ^. The most interesting are the diagrams of 
Fig. Pd), c) which are of the type of renomalons and should behave to k-th order of 
perturbation theory as k\. 

In order to calculate D'^{n;p) for the unbroken 0{N) symmetry at large A^, we first 
derive the proper Schwinger-Dyson equation substituting F = ^''(0, —ir') in Eq. ( |2.10| ) 
and symmetrizing according to Eq. ( |4.17| ). The resulting equation can be written as 




ir)|0) =-S{t-t'). (4.18) 



conn 



This equation involves a new amplitude 

D\r-p,,p,) = ^j ^^e^^-+^^^^-(-+-)^(H|0''(x)/(y)0"(O,-.r)|o)_ (4.19) 

with three external lines having nonvanishing spatial momenta (which is a generating 
function for D"'{n;pi,p2)), where a is the 0{N) index of an outgoing particle with the 
spatial momentum pi + p2 while the averaging over the 0{N) indices of two incoming 
particles with the spatial momenta pi and p2 is performed. 
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Figure 6: The graphic representation of Eq. ( 4.20| ). The vertex on the l.h.s. represents 
D"'{t;pi,P2), the circles with two hnes represent D'^{t;p) or D^{t;p) and the 
tadpole represents <I>'^(t). 



The Scwhinger-Dyson equation which is satisfied by D"'{t;pi,P2) is now of the vertex 
type rather than of the propagator type as before. It is depicted graphically in Fig. ^ 
and reads analytically as 

^-D'{t;p,+P2)+ j ^^D\r-k,p,+p2-k) 



X 



(4.20) 



Its solution for D"'{T;pi,p2) is given by 

XND^ir; p,)D^{r; j92)<l>"(r)Z}^(r; p, + p^) 



D''{t;pi,P2) 



1-XNJ (0/^^(r; k)DTiT; p,+P2- k) 



(4.21) 



The propagator D'^{t;p) is determined by Eq. ( p.l8| ) and reads explicitly 



D^{T;p) = J^e<^'-^)GUr,r') 



2i 



^n=0 V 8 



-2mT 



+ 



2i 



1 



1 Ml „ -2m7 
1 8 



€m—2nm^ 
{€—2nm)^—uj^ 

n 



-1 I 2m(m+e) 



(4.22) 



where the integral over dr is understood as an analytic continuation from imaginary r 



(cf. (|2.8|) ). The expression ( [4.22| ) can be simplified using the formula 



oo 



E 



u ^ 2Fi(1,m; 1 + u;z). 



The calculation of the renormalon contribution would consist of three steps. 
1. Calculate the integral / ^^D'^{t; k)D'^{T;p- k). 



(4.23) 



2. Substutute the (proper Fourier transformed) expression (4.21) into Eq. ( [4.18 ) and 
solve it for the diagonal resolvent G^{t,t). 

3. Integrate dcu to ^uP- — w? G^{t, t). 



II This equation differs at the tree level from the one ofRef. Il) since we use exact propagators which 
are "dressed" by emitting on-mass-shell particles. 
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a) b) c) 

Figure 7: The tree diagrams for a 2— >3 amplitude in matrix (f)^ theory at large A^. The 
continuous lines are associated with propagation of matrix indices. The contri- 
butions of the diagrams a), b) and c) equal, respectively, 1/4, —1/6 and 1/24 for 
the given on-mass-shell kinematics when the momenta of the incoming particles 
(in the units of mass) are (3/2, a/5/2, 0, 0) and (3/2, -\/5/2, 0, 0), while that of 
the each produced particle is (1,0,0,0). The diagrams b) and c) have an extra 
combinatorial factor 2 each. The sum of three diagrams vanishes which illustrates 
the nullification of 2— >3 on mass shell at the tree level. 

4.3 Matrix Higgs field 

Another interesting model to investigate is the case of the matrix Higgs field which is 
described by the Lagrangian 

£ = 1 tr {d,<pf - !^ tr 02 - ^ tr 03 - ^ tr 0\ (4.24) 

where (j)^^{x) is generically NxN Hermitean matrix. The model greatly simplifies as 
N ^ oo at fixed X^N or A4A^ when only the planar diagrams survive similar to the 
't Hooft large- limit of QCD. 

The nullification of on-mass-shell amplitudes holds in the case of the matrix cubic 
interaction for 2— >-3 at large at the tree level. The cancellation of diagrams is illustrated 
by Fig. 1^. This looks quite similar to the 0(A^)-vector case. 

The large- limit of the matrix 0^ theory while being simpler than the A^ = 1 case 
to all loops is, however, quite nontrivial like the large-A^ limit of QCD. Some of the 
rescattering diagrams are now included in the large- A^ amplitudes so they should not look 
like the tree-level ones. The crucial problem is to find a simple set of observables like 
tr0'^ for ordinary matrix models to close the Schwinger-Dyson equations for the given 
kinematics at large A^. 
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5 



Conclusion and some further problems 



• Quantum field theory is not yet completely investigated for large multiplicities of 
identical particles (n ~ 1/A). 

• Both conventional perturbation theory and the expansion break down at n ~ 
1/A or ~ A^. 

• While for n~ 1/A or n~iV the problem becomes semiclassical, there are subtleties 
in application of semiclassical technique. Alternative methods could be useful. 

• The mechanism of unitarity restoration is not yet completely work out (while there 
is no problems with unitarity at large N). An exponentiation of the corrections in 
n^A or /N looks promising. 

• Cross-sections of multiparticle production may become large at very high energies 
E ~ m/A while smaller than the unitary bound. 

• The nullification survive at large A^ to all loops (factorization is crucial for the 
reduction to a quantum mechanical problem at the tree and one-loop levels or at 
large A^). 

• Which dynamical symmetry is behind the nullification is not yet found out and 
higher conservation laws which restrict multiparticle production are not yet con- 
structed. 

• Drastic simplifications occur for the given kinematics (in contrast to the case when 
all momenta are off mass shell). The threshold amplitudes are described by an 
integrable lower dimensional problem. 

• One might try to solve more complicated D = A models imposing this kinematics. 

• The developed technique might be applicable to some other problems. 

• The nullification could have some phenomenological consequences for the standard 
model restricting the ratio of masses. 

• An extension to nonvanishing spatial momenta would be very interesting, in partic- 
ular, for matching |Q the threshold behavior of multiparticle amplitudes with the 
standard multiperipheral picture of high energy collisions. 
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